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摘要 : 空间 探测 器 对 X 射线 天 体 源 的 有 效 观 测 时 间 主 要 受到 空间 环境 因素 的 制约 。 制 约 有 
效 观 测 时 间 的 主要 因素 包括 太阳 避免 角 ， 地 球 对 源 的 遮挡 及 南大 西洋 异常 区 等 ; 然而 
卫星 处 于 一 些 高 粒子 本 底 区 域 ， 太 阳 与 地 球 之 间 的 光照 区 及 视 场 指向 接近 地 球 时 ， 本 
底 水 平 很 高 而 且 难 以 确定 ， 这 些 时 间 的 数据 也 是 难以 使 用 的 。 本 文 利 用 orbitTools xf 
数 库 预测 轨道 ， 利 用 HEAsoft 下 的 attitude 函数 库 计 算 空 间 环境 变量 ， 利 用 这 些 空间 
环境 变量 对 源 的 观测 时 间 进 行 估计 ， 并 通过 与 实际 观测 比较 ， 证 明 本 方法 估计 的 观测 
时 间 能 够 比较 好 的 与 实际 情况 一 致 。 
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Abstract: The observation constraints for astronomical satellite are always due 
to space environment, such as Sun angle constraint, the time occulted by the 
Earth and the time passing through South Atlantic Anomaly. There are also some 
time periods, such as the time in high background region, the time when 
satellite between the Sun and the Earth, and the time when field of view of the 
satellite is close to the surface of the Earth, which can' t be used for 
scientific analysis since the background level is very high or the background 
model has a great of uncertainties. In this work, we establish an interface for 
these variations to estimate observable time of X-ray sources, and use 
orbitTools library to generate orbital data and HEAsoft' s attitude library to 
calculate these variations. Compared with the actual observation, the predicted 
time of this work is consistent well with the time of observed data. 
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0 引言 

便 X 射 线 调制 望远镜 (Hard X-ray Modulation Telescope, HXMT) 卫星 DPI 在 2017 年 6 
月 发 射 升 空 ， 开 启 了 中 国 在 六 射 线 波段 对 天 体 源 的 观测 研究 。 目 前 爱 因 斯 坦 探 针 卫 星 外 
(EP) 己 经 正式 并 项 ， 并 希望 于 2022 EFFE, MIRE X 射线 时 变 与 偏振 天 文 台 (eXTP) 
也 在 预 研 阶段 。 这 些 卫 星 的 轨道 高 度 大 多 处 于 地 球 内 辐射 带 范 围 息 ， 卫 星 载荷 容易 受到 南 
大 西洋 异常 区 “SAA) 及 地 磁场 的 影响 外 ,使 得 卫星 获取 的 数据 质量 很 差 ( 本 底 高 且 形 状 
难以 确定 ) ， Rc 区 数据 ， 比 如 ， 多 数 卫 星 进 入 SAA 区 后 停止 采集 数据 或 者 直接 关 
机 。 同 时 这 些 卫星 的 一 个 重要 任务 是 对 源 进行 定点 观测 ， 而 源 的 可 见 时 间 主 要 受到 太阳 及 
地 球 的 影响 : 卫星 的 热 控 不 允许 载荷 指向 太阳 而 是 需要 避 开 太阳 一 定 角度 ， 比 如 HXMT 
有 遮阳 板 ， 并 同时 要 求 卫星 指向 避 开 太阳 ， 这 决定 着 观测 源 需 要 远离 太阳 位 置 ， 地 球 会 距 
挡 卫 星 视 场 ， 使 源 的 光子 无 法 被 载荷 获取 。 

过 去 ， 天 体 源 观测 时 间 的 估计 主要 集中 在 对 源 可 见 时间 的 估计 钼 上， 通常 的 做 法 是 引 
入 观测 约束 ， 可 见 时 间 即 为 经 过 约束 条 件 筛选 之 后 的 时 间 ， 这 些 观 测 约束 包括 地 球 对 视 场 
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的 遮挡 ， 太 阳 和 月 球 对 视 场 的 遮挡 以 及 对 热 控 的 约束 ， 和 经 过 南大 西洋 异常 区 的 时 段 。 实 
际 上 这 些 约束 是 不 够 的 ， 源 的 观测 时 间 还 要 排除 一 些 特定 的 区 域 ， 比 如 高 粒子 本 底 区 域 ， 
而 这 些 约束 一 般 在 数据 分 析 时 才 考 虑 ， 比 如 对 低 轨 天 文 卫星 Swit, Suzaku, EENE 
布 的 数据 分 析 手 册 中 均 给 出 地 磁 截 止 刚度 CCOR) 外 的 筛选 范围 《以 排除 高 本 底 区 域 ) 。 
这 样 ， 一 个 X 射线 天 体 源 的 观测 时 间 佑 计 需 要 从 数据 质量 〈 包 括 本 底 水 平 及 本 底 形 状 不 确 
定性 ) 及 可 见 时 间 两 方面 着 手 。 但 是 数据 质量 的 约束 条 件 与 载荷 的 工作 机 制 及 所 在 能 区 是 
有 关系 的 ， 不 同 的 载荷 需要 不 同 的 约束 条 件 ， 比 如，HXMT 的 低能 望远镜 中 ， 当 卫星 处 于 
太阳 和 地 球 之 间 时 ， 其 SCD (Swept Charge Device) 阵列 会 受到 不 同 程度 的 漏 光 ， 所 以 需 
要 排除 此 部 分 时 间 。 

本 文 主要 关注 低 轨 天 文 卫 星 所 在 的 空间 环境 对 源 观测 时 间 的 影响 ， 并 从 可 见 时 间 和 数 
据 质量 两 个 方面 给 出 估计 方法 ， 这 将 有 助 于 观测 计划 的 制定 ， 提 高 天 文 卫星 的 利用 ， 促 进 
科学 成 果 的 获取 。 
1 天 文 卫 星 影 响 因 素 分 析 

观测 时 间 估 计 中 ， 我 们 只 考虑 空间 环境 因素 的 影响 ， 实 际 上 卫星 指向 的 晃动 及 载荷 温 
度 的 变化 也 会 带 来 观测 时 间 的 下 降 ， 但 这 些 因素 与 望远镜 结构 性 质 有 关 ， 不 在 本 文 讨论 范 
围 内 。 

1.1 地 球 对 观测 时 间 的 影响 

这 主要 表现 在 三 个 方面 : 〈1) 地 球 将 视 场 遮挡 ， 源 光子 不 能 进入 视 场 ，(2) 卫星 视 
场 靠 近 大 气 层 时 ， 散 射 的 太阳 X 射 线 光 子 或 者 可 见 光 进 入 视 场 内 ， 导 致 本 底 难 以 估计 ; 
G) 源 的 光子 经 过 大 气 层 吸收 也 会 导致 辐射 的 畸变 。 

一 般 用 ELV，DYE_ELV，NTE_ELV 表示 地 球 的 影响 ， 它 们 分 别 为 视 场 指向 与 地 球 最 
小 的 角度 ， 指 向 与 地 球 明亮 区 〈 太 阳 照 射 的 一 面 ) 最 小 的 角度 ， 和 指向 与 地 球 阴 暗 区 最 小 
的 角度 ， 如 果 后 两 者 不 存在 ， 则 一 般 设 为 120° 或 者 200"。ELV Æ DYE ELV 一 般 越 大 ， 地 
球 对 观测 的 影响 越 小 ， 但 是 观测 时 间 也 会 降低 。 对 于 和 射线 探测 器 ， 如 果 探 测 能 区 比较 低 ， 
ELV 可 以 适当 大 一 些 以 降低 对 低能 端 能 谱 的 影响 ， 时 变 分 析 可 以 要 求 ELV 适当 降低 ; 
DYE ELV 主要 反映 太阳 弥散 和 射线 本 底 和 可 见 光 对 仪器 的 影响 ， 这 对 CCD 型 探测 器 影响 
比较 大 。 

1.2 太阳 及 月 球 的 影响 

太阳 对 卫星 的 影响 是 巨大 的 ， 主 要 影响 热 控 ， 而 且 太 阳 X 射线 也 会 对 和 射线 探测 器 有 
影响 ， 尤 其 是 对 低能 区 探测 器 影响 比较 大 。 一 般 约束 卫星 的 姿态 使 其 上 面 的 载荷 偏离 太阳 
视 场 指向 与 太阳 中 心 的 夹 角 (SUN_ANG 表示 ) 越 大 越 好 。 
HER (MOON ANG， 指 向 与 月 球 夹 角 表 示 ) 会 散射 太阳 X 射 线 和 可 见 光 ， 对 视 场 有 


c 0.5" 立 体 角 的 遮挡 影响 ， 一 般 越 大 越 好 。 
O 1.3 南大 西洋 异常 区 的 影响 


SAA 区 有 很 高 的 粒子 本 底 背 景 ， 对 载荷 的 损坏 很 大 ， 一 般 载 荷 经 过 它 时 要 关机 。SAA 
的 边界 不 明显 ， 接 近 SAA 时 ， 粒 子 本 底 会 逐渐 变 大 ， 这 一 段 时 间 也 是 不 可 用 的 ;离开 
SAA 区 时 ， 卫 星 的 一 些 材料 因为 受到 SAA 区 粒子 的 影响 而 活化 加 重 ， 使 得 活化 本 底 很 高 ， 
然而 活化 衰减 也 比较 快 ， 从 出 SAA 到 活化 衰减 到 一 定 程 度 所 用 的 时 间 往 往 也 是 不 可 用 的 。 

SAA FLAG, T SAA 和 TN_SAA 色 表示 卫星 是 否 在 SAA， 出 SAA 的 时 间 长 度 ， 和 进 
下 一 次 SAA 的 时 间 长 度 ， 一 般 此 三 量 表示 SAA 的 影响 ， 需 要 注意 这 两 个 量 与 指定 的 SAA 
区 大 小 有 关 。 

1.4 高 粒子 本 底 区 域 的 影响 

除了 SAA 区 的 影响 ， 一 些 高 粒子 本 底 区 域 ， 比 如 高 纬度 地 区 ， 对 数据 质量 影响 也 很 
X. COR 是 描述 地 球 磁场 阻挡 宇宙 线 能 力 的 量 ， 粒 子 要 进入 地 磁场 中 的 某 一 位 置 就 必须 大 
于 此 处 的 地 磁 刚 度 ， 它 一 般 用 最 小 动量 表示 。 因 为 宇宙 线 会 引起 仪器 本 底 增多 ， 所 以 COR 
比较 低 的 地 方 ， 本 底 会 比较 高 。 另 一 方面 ， 一 般 不 能 直接 测量 载荷 的 本 底 ， 而 是 利用 模型 
估算 本 底 ，COR 比较 低 的 地 方 ， 因 受 低 能 电子 、 质 子 影 响 大 ， 本 底 形 状 不 确定 性 也 就 很 大 。 
星 处 于 高 的 粒子 本 底 区 域 ， 光 照 区 及 视 场 接近 地 球 边 缘 时 ， 本 底 水 平一 般 比 较 高 ， 
且 形 状 难以 估计 ， 使 得 卫星 采集 的 数据 难以 使 用 ， 这 会 使 得 观测 时 间 的 估计 与 真实 差别 比 
较 大 。 
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2 观测 时 间 估 计 方 法 

我 们 的 目的 是 估计 某 一 源 或 者 多 个 源 在 某 段 时 间 内 的 观测 时 间 。 这 一 实现 将 有 助 于 对 
多 个 观测 源 进 行 观 测 效率 〈 观 测 时 间 与 总 时 间 之 比 ) 优化 ， 从 而 制定 高 效率 的 运行 计划 。 

我 们 利用 orbitTools00 软 件 包 对 卫星 轨道 进行 预测 ， 利 用 HEAsoft"?' F IJ attitude 软件 包 
对 空间 环境 进行 参量 估计 。orbitTools 软件 包 基 于 SGPA/SDPA 模型 ， 分 别 对 低 轨 和 高 轨 卫 
星 进行 轨道 预测 ， 输 入 为 两 行 轨 道 根 数 及 参考 时 间 。 两 个 软件 包 均 采用 cc++， 方 便 集 成 。 

2.1 源 信息 及 约束 条 件 的 定义 

观测 时 间 估 计 中 源 信息 其 实 不 是 特别 重要 ， 重 要 的 是 卫星 指向 信息 ， 这 在 偏 轴 观 测 中 
尤为 特别 ， 所 以 这 里 指 的 源 信 息 其 实 是 卫星 指向 信息 。 我 们 定义 在 J2000 F: z 轴 为 卫星 指 
向 方向 ，y 轴 为 卫星 某 一 轴 〈 比 如 卫星 本 体 坐 标 系 下 x 轴 ) 指向 太阳 的 方向 。 

为 了 使 处 理 流 程 是 按时 序 进行 的 ， 我 们 将 SAA 区 域 扩 大 ， 从 而 简化 T SAA 和 
TN_SAA。 实 际 中 单纯 通过 硬件 触发 识别 SAA 是 存在 问题 的 ， 卫 星 靠近 SAA 区 域 (而 不 
AT SAA) ， 载 和 荷 计 数 也 会 发 生 畸 变 ， 最 好 的 方法 是 积累 大 量 数据 之 后 通过 分 析 数 据 进 而 
定义 SAA。 

这 样 源 及 约 R 信 息 为 : 
RA, DEC, ELV, DYE ELV, COR, SAA, MOON ANGLE, SUN ANGLE, ， 其 中 
RA, DEC 为 卫星 指向 位 置 ，SAA 表示 几 种 SAA 的 模型 之 一 ， 模 型 可 以 由 用 户 定 义 。 


22 轨道 及 运行 参数 定义 
我 们 利用 orbitTools 软件 包 进行 轨道 预测 ， 它 的 输入 为 两 行 轨道 根 数 《TLE》。 
运 f E 数 包 in 


tl, t2, stepsec, year, month, day, hour, min, second, mjd0, output. $F, tl, (X 
示 自 参考 时 刻 秒 的 累计 量 ,， 估计 的 观测 时 段 从 开始 至 也 结束 ; 参考 时 刻 由 
year, month, day, hour, min, second 表示 ， 比 如 对 于 HXMT， 它 们 分 别 为 
2012，1，1，0，0，0，second 可 以 为 浮 点 数 ; mjd0 为 参考 时 刻 的 MJD 时 间 ， 对 于 HXMT 
为 55927; output 为 输出 文件 的 名 称 ， 本 设计 中 它 的 格式 为 天 文 数据 格式 FITS; stepsec 
为 时 间 步 长 (整数 ) ， 单 位 为 秒 ， 表 示 每 积累 stepsec 秒 计算 一 次 观测 效率 。 
轨道 预测 中 所 有 时 间 均 为 相对 于 参考 时 刻 的 秒 累计 量 《UTC 时间) ， 但 是 在 计算 ELV 
等 量 的 时 候 需 要 将 这 些 UTC 转换 为 MJD 时间， 需要 考虑 的 是 两 者 之 间 存 在 跳 秒 。 
基于 目前 的 一 些 商业 软件 〈 比 如 STKO 或 者 一 些 机 构 ( 中 国 科学 院 空间 中 心 ) 专门 可 
以 对 轨道 进行 预报 ， 轨 道 也 可 以 使 用 第 三 方 ， 这 样 程序 最 好 也 可 以 读 入 FITS 格式 的 轨道 文 
件 ， 其 需要 包括 时 间 CTIMEO MME (X, Y, Z2) 信息 ， 而 且 需 要 在 J2000 坐标 系 下 ， 这 
与 orbitTools 产生 的 轨道 坐标 系 一 臻 《当然 ， 也 可 以 根据 时 间 将 J]2000 转换 到 地 固 坐 标 系 ) 。 

2.3 实现 策略 

软件 的 参数 读 取 采用 文件 读 取 方 法 ， 将 源 及 约束 条 件 ， 两 行 轨道 根 数 及 运行 参数 分 别 
设计 成 配置 文件 〈 源 及 约束 条 件 文件 ，TLE 文件 和 运行 参数 文件 ) 。 如 果 要 估计 多 个 源 同 
时 段 的 观测 时 间 ， 可 以 将 多 个 源 位 置 及 其 各 自 约 束 写 入 到 源 及 约束 条 件 文件 中 。 如 果 轨 道 
预测 不 采用 orbitTools 库 ， 则 直接 读 入 轨道 数据 ， 对 于 轨道 数据 点 之 外 的 点 用 线性 插值 方法 
计算 得 到 其 位 置 。 
程序 首先 初始 化 ， 读 入 上 述 配 置 文件 ， 找 到 UTC 与 MJD 对 应 关系 ， 指 定 效率 计算 的 
起 止 时 间 ， 利 用 TLE 及 MIDO 初始 化 orbitTools 各 个 模块 ， 同 时 设置 约束 条 件 的 取 值 ， 然 
后 计算 一 次 时 间 步 长 (stepsec〉 内 的 轨道 点 及 其 时 间 〔 即 1s 一 个 轨道 点 ) ， 计 算 这 些 点 的 
ELV, COR 等 值 并 与 约束 条 件 比 较 ， 统 计 符合 约束 的 点 数目 。 时 间 步 长 中 的 所 有 点 判断 完 
成 后 ， 观 测 效率 即 为 符合 约束 的 点 与 总 的 点 〈 即 stepsec) 比值 。 时 间 步 长 如 果 比 较 小 〈 比 
如 1s) ， 则 计算 的 观测 效率 要 么 为 0 要 么 100%。 时 间 步 长 应 该 从 fl 开始 ， 直 到 超过 人世 为 
IE; 最 后 将 每 段 观 测 效 率 及 其 时 间 (中 间 时 间 表 示 〉 写 入 文件 。 

COR 与 卫星 轨道 有 关 ， 我 们 采用 了 IGRF2005 外 产生 的 结果 ， 该 结果 的 纬度 范围 比较 广 ， 
这 对 一 些 倾角 比较 大 的 卫星 HXMT (卫星 倾角 43") 比较 合适 。 
3 预测 结果 及 讨论 

我 们 采用 HXMT 卫星 的 高 能 望远镜 在 2017 年 8 月 底 的 一 次 Crab 观测 结果 与 本 方法 进 
行 比较 ， 同 时 要 求 ELV>10"，COR>8GV，RA，DEC 分 别 为 255.705667 和 -48.7896。 1 
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为 实际 结果 与 预测 结果 的 比较 ， 横 坐标 为 时 间 〈 单 位 为 秒 ) ， 自 上 而 下 分 别 为 Crab 原始 的 
观测 光 变 Nal 事例 ， 扣 除了 地 球 遮挡 及 SAA DO ， 利 用 HXMT HIP Uti y Br ATE UG 
后 的 光 变 ， 本 方法 预测 的 观测 时 间 (观测 效 率 ) 。 原 始 光 变 来 自 HXMT 提供 的 1 级 数据 产 
品 ， 经 过 筛选 后 有 相当 多 的 时 间 被 扣除 ， 而 本 方法 计算 的 效率 或 者 说 可 用 时 间 与 筛选 后 的 
结果 比较 一 致 。 需 要 说 明 的 是 本 方法 的 姿态 是 按 卫 星 本体 坐 标 系 下 x 轴 指向 太阳 ， 并 随 着 
时 间 姿 态 会 有 所 变化 ， 而 实际 观测 中 ， 姿 态 一 旦 确定 ， 整 个 观测 是 不 变 的 。 
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图 l, crab 观测 结果 与 预测 效率 的 比较 。 自 上 而 下 分 别 为 Crab 原始 观测 结果 ， 经 过 数据 得 
选 后 的 结果 及 预测 的 效率 。 
Fig1，Comparison between Crab observation data and the predicated efficiency. The top panel 
shows the Crab Nal light curve, the middle panel shows light curve after event screening by HXMT 
user analysis software, while the bottom panel shows the predicated efficiency by this work. 


可 见 时 间 的 估计 不 同 于 观测 时 间 ， 观 测 时 间 还 需要 关注 数据 质量 ， 特 别 是 对 于 一 些 低 
能 探测 器 ， 比 如 CCD 型 探测 器 ， 它 可 能 存在 漏 光 ， 这 使 得 DYE_ELV 角度 非常 重要 ， 而 且 
还 考虑 了 COR 等 的 影响 ， 使 得 效率 直接 反映 出 数据 的 使 用 情况 。 

但 是 轨道 预测 的 精度 受到 TLE 的 影响 ， 但 另 一 方面 即使 TLE 准确 ， 卫 星 因为 受到 大 气 
影响 ， 地 球 重力 不 均匀 等 因素 ， 长 时 间 预 测 也 是 不 准 的 。 
目前 卫星 能 源 供给 比较 充裕 ， 一 般 情况 下 可 以 允许 卫星 在 一 轨 中 多 次 调 向 ， 这 对 提高 
卫星 的 观测 效率 是 相关 重要 的 ， 以 本 方法 为 基础 ， 在 一 轨 中 考虑 如 何 提高 观测 效率 ， 特 别 
是 对 一 些 机 动 性 能 比较 好 的 小 卫星 ， 比 如 EP， 是 至 关 重 要 的 。 

本 文 没有 考虑 太阳 观测 卫星 的 观测 时 间 估 计 ， 事 实 上 太阳 也 是 一 颗 和 射线 天 体 源 。 另 
外 本 文 也 没有 考虑 地 面 大 功率 设施 对 卫星 数据 的 影响 。 
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